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We discuss the effective field theory approach to nuclear phenomena with typical 
momentum below the pion mass. We particularly focus on the three body problem. 
In the J = 3/2 channel, effective field theory is extremelly predictive and we are 
able to describe the nucleon-deuteron scattering amplitude at the percent level 
using only two parameters obtained from nucleon-nucleon scattering. We briefly 
comment on the J = 1/2 channel and the issues related to the non-perturbative 
renormalization of the three body force. 



1 Introduction 

Most of the early work on the use of effective field theory in nuclear physics 
concentrated on the momentum scale p ~ m^. There is a good reason for 
that: most of nuclear phenomenology is at this scale (the Fermi momentum in 
nuclear matter, for instance, is at the order of 2m7r). This talk, however, will 
concentrate on the scale p ~ 1/a (a is the nucleon-nucleon scattering 

lenght), that will be refered to as the " very low " momentum scale. This 
scale region is interesting for three reasons. The first reason is that this scale, 
even though not typical in nuclei, is relevant for a corner of nuclear physics. 
This comes about because the deuteron is an anomalously shallow bound state, 
with a binding energy of the order of \/Ma^ , what sets the momentum scale of 
deuteron physics to be p 1/a ^ m^. A second reason is that, since pions can 
be integrated out and don't appear as explicit degrees of freedom, a number of 
issues related to regularization and renormalization, many of them extensively 
discussed in this workshop, can be addresed in a simpler situation. The third 
motivation to consider the " very low " momentum scale is related to recent 
work on the " low " momentum scale p ^ to^ . At this higher momentum scale 
pions have to be included explicitly and in previous power counting schemes the 
whole ladder of pion exchanges contributed at leading order to nucleon-nucleon 
scattering. This is a source of tremendous complication, both conceptual and 
practical. Recently a new power counting scheme was suggested pp in which 
pion exchanges are perturbative and, in particular, the leading order contains 
no pion exchanges. Thus, all results obtained discussed here using the effective 
theory appropriate to the p <C m-^ scale will automatically equal the leading 
order result obtained from this new power counting of valid at the p ~ TO,r 
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scale. 

In the absence of explicit pions, the effective theory method applied to 
nucleon-nucleon scattering simply reproduces the effective range expansion. 
This is phenomenologically correct, but hardly interesting. The first nontrivial 
application is the three body problem, that is the main subject of this talk. 



2 The very low energy effective theory 

For momenta much smaller than the pion mass, the only degrees of freedom 
that need to be included explicitly are the nucleons. The effect of pion ex- 
changes, A's, heavier meson, etc., are all implicitly included in the coefficient 
of local operators in the lagrangean. We will concentrate later on the three 
body problem in the J = 3/2 channel. In this channel all nucleon spins are 
parallel thus all s-wave two nucleon interactions are in the spin triplet channel. 
Restricting ourselves to spin triplet interactions, the most general lagrangean 
is : 



+C2 [(A^V2CTCT2VAr)(7VV2CTCT2V7V) - 3(7V V2CTCT2iV) (TV ^2(7(72 V^iV) + h.c] 
+ ..., (1) 

where AI is the nucleon mass, C„ are constants related to the two-body force 
terms containing n derivatives, and the dots stand for higher-order terms in- 
cluding relativistic corrections, higher-derivative terms, three-body forces, etc. 
Terms describing the spin singlet interactions have a similar form but will not 
be needed here. The constants C„ are determined by nuclecm-nucleon scatter- 
ing data. As it was explained in more detail in van Kolck'sEI talk it turns out 
that, using dimensional regularization and minimal subtraction, Cq ^ a/M, 
C2 ~ ?'o(?'oa)/A^i Ci ^ ro{roa)'^ /M + . . . and so on (ro is the effective range 
and ellipses stand for terms suppressed by powers of rg/a ). The leading pieces 
in each one of these terms form a geometric series that can be convenieptly 
summed to all orders by the introduction of a field of baryon- number two a 



-|((? •iVaa2T27V + h.c.) + ... (2) 

If the dibaryon field d is integrated out, the Lagrangian (|l|) is recovered as long 
as A and g are appropriate functions of Cq and C2. This resummation is by 
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no means necessary, since for momenta of the order p ^ 1/a the resummcd 
terms are subleading, but it is just a convenient way of computing higher-order 
corrections. 

The numerical values of g and A can be determined if we consider the 
dressed dibaryon propagator (Fig. |l|). 





Figure 1: Dressed dibaryon propagator. 

The linearly divergent loop integral is set to zero in dimensional regularization 
and the result is 



'^^P) = ; = • (3) 

4RI '-^'2tt 



This propagator is, up to a constant, the scattering matrix of two nucleons in 
the ■^S'l channel, 

47r 1 

where jM is the energy in the center-of-mass frame. This result is just the 
familiar effective range expansion, from what we can infer thejproper values 
for the constants g and A. Using a — 5.42 fm and rg = 1.75 fmiD, we find 

5' = 17% = 1-6 • 10-' MeV-\ (5) 
A=M!7;7-8.7MeV. (6) 

From Eqs. (|), (|), and M) we see why it is necessary to resum the bubble 
graphs in Fig. ^ to all orders for p ^ 1/a: the term in the square root coming 
from the unitarity cut is of the same order as A. On the other hand, as 
mentioned before, the kinetic term of the dibaryon is smaller than the other 
terms in (^) and is resummed for convenience only. Notice that the propagator 
has two poles, one at ~ p'^ /4M — B (the deuteron pole), another at 
= p"^ /AM — Bdeep (unphysical deep pole), and a cut along the positive real 
axis starting at p° = p'^/AM. 
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3 Power counting for neutron-deuteron scattering 

Let us now turn to neutron-deuteron scattering. We want to identify which 
graphs give the leading contributions in an expansion on powers of the typical 
momentum p 1/a over the scale of the physics not explicitly included in the 
effective theory m^r ^ 1/t'o- 

The simplest diagram contributing to neutron-deuteron scattering is the 
first diagram on the right hand side of Fig. 0. 




Figure 2: Dressed dibaryon propagator. 

For momenta of the order of p ~ 1/a it gives a contribution of the or- 
der of (7^/(p^/Af) ~ a? /Mtq. The one- loop graph mixes different orders of 
the expansion, since it involves the dibaryon propagator ^^/(A /M) ~ 
(a/M)(l -I- ©(ro/a) -f . . .). It also involves a loop. The contribution com- 
ing from loops can be estimated by rescaling all internal energies by p^/M, 
, where p is some typical external momentum. The result is M^p"^ , where 
n =number of nucleons propagators minus the number of loops and m can be 
determined by dimensional analysis. Using this rule we can easily estimate the 
one loop graph in Fig. |to be g'^g'^ / {I^+p^ /M) ^ {a'^ /Mro){l + 0{ro/a) + . . .). 
Similarly we can see that the remaining graphs in Fig. |^ give contributions of 
the same order, which means that an infinite number of diagrams contribute 
to the leading orders. 

Other contributions are suppressed by at least three powers of rg/a or 
pro- For instance, the effect of the subleading (not resummed) piece of C4 is 
to generate the shape parameter fc'^) term in the effective range expansion 
of the nucleon-nucleon interaction. Its typical size is ~ k^rQ compared to the 
leading piece 1/a and is thus also suppressed by (rg/a)^. Likewise, p-wave 
interactions, unaffected by the existence of a shallow s-wave bound state, arise 
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from a term in the lagrangean with two derivatives and a coefficient of the order 
^ 1/MmJ. We conclude then that a diagram made out of the substitution of 
one of the dibaryon propagators in a diagram in Fig. |^ by a p-wave interaction 
vertex would be suppressed by (ro/a)"^ in comparison to the leading order. 
Three-body force terms have to contain at least two derivatives since in the 
J = 3/2 channel all the spins are up and Fermi statistics forbids the placement 
of all three nucleon in a s-wave. The natural size of the coefficient of the six 
nucleon, two derivative term is 1/Mml. A dia gram including this three body 
force includes two nucleon loops, each one giving a Mp factor, two dibaryon- 
two nucleons vertices and a three body force vertex, for a total contribution 
of the order oi p'^ /Mm^ro ^ a'^/Mro{ro/a)^. Thus contributions coming from 
the three-body force are suppressed in relation to the leading order graphs by 
(r-o/a)6. 



4 Summing the leading graphs 

From the results of the previous section we see that a calculation accurate up 
to corrections of order (rp/a)'^ is possible by summing the diagrams of Fig. ^ 
Fortunately, the interaction mediated by the s-channel dibaryon generates a 
very simple, local and separable potential between nucleons. It is well known 
that the three-body problem with separable two-body interactions reduces to 
an equivalent two-body problem. In our case the equation to be solved can be 
read off Fig. 0, and an integration over the energy inside the loop givesQ 



i{p, fc) = ^ 

{p~k/2Y + MB 



dH i{l,k) 



1 

p^f-p + p^ - IP + mb' 

where B is the deuteron binding energy and i is a correlator with two dibaryon 
and two nucleons external legs represented by the blob in |^. One could have 
computed the deuteron-nucleon amplitude from the correlator involving any 
other operator with the deuteron quantum numbers instead of the dibaryon 
used here. The results for the (on shell) deuteron-nucleon scattering are ob- 
viously independent of this choice as long as the proper wave function renor- 
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malization is performed. The scattering matrix is given by 

t{p) = Vz i{p,p) Vz, 

where Z is defined by 

1 



Z = iSo(p)\po\e/iSip)\po\e 



1 - 



(8) 



(9) 



4-kVMB 

It is convenient then to define the off'-shell scattering ampHtude 



t{p, k) 



satisfying 



'-^im^ + hNW - k') +MB- VMB) 



t{p,k), 



(10) 



dH 1 



t{p, k) 



{p- k/2y 



MB 



^2 _ ^2 _ jg 

t{l,k) 



(11) 



(27r)3 12 ^i.^j^p_^k'^j^MBp-k'^-ie 



that reduces to the scattering ampHtude t{p) on shell 

t{p) =t{p,p). 



(12) 



Since we are interested only in s-wave scattering, we should project this 
equation into its L = component. The result is 



1 



a{x,y) = - —In 
xy 



TTX 



J dz zln ^ 



2;2 + + 1 - |y2 _^ xz 



z^+x^ + l-h^- 



ix + y/2)^ + l 
{x - y/2f + 1 

a{z,y) 



XZ J z"^ - y"^ - le 



(13) 



where wc use the dimcnsionless quantities x — p/yMB, y = k/yMB, z 



1/VMB, and a{x,y) = ^^^tL=Qip,k), and 77 = y/MBro/2. For finite values 
of k this equation is complex even below threshold (3fe^/4 = B) due to the 
ie prescription. The numerical solution of the equation above is trivial in 
terms of computer power. The only subtle point is how to deal with the ie 
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prescription that appears in eq.(^3|). One way of dealing with that is to use 
the real X-matrix defined by 



i{x,y) 



which satisfies the equation 



2 -p 

K{x, y) = -h{x, y,y) / dz z h{x, y, z)^ :;K{z, y), 

TT Jo z ~ y- 



(14) 



(15) 



with 



h{x,y, z) = 



-In 



+ + 1 - ly^ + xz 



xzf{x,y) \z'^ + x"^ + 1 - jy"^ - xz J ' 



f{x,y) = 2 



^|(x2-y2) + i + i_ 
The phase shifts can be obtained directly form the on-shell iiT-niatrix 



(16) 



fccot(5 



y/MB 



ITS' Vmb 



(17) 



Defining f{x,y) by the equation 
h{x,y,y) 2 



h{y,y,y) tt 
the on-shell if-matrix can be obtained by 



dzz [h{x,y,z)~— -h{y,y,z} 

Hy,y,y) 



fjz.y) 

9 9 ' 



K[y,y) ^ -h{y,y,y) 



2 

1 + - 

TT 



dz {z^h{y, y, z)f{z, y) - y'^h{y, y, y)f{y, y)) 



9 9 



(18) 



(19) 

-1 



Rewriting Eq. ( [l^ ) this way greatly simplifies its numerical solution, for 
now the integrand is regular and the principal value can be dropped from Eqs. 
(||)and(|l|). n ^ 

We have solved Q Eqs. ( [l8| ) and ([l^ ) numerically and the result for the 
phase shifts for energies up to the break-up point is shown in Fig. ^.j-The data 
points at finite energy were taken from the phasej-|Shift analysis inl3 and the 



much more precise (nearly) zero-energy point fromB. Also plotted is the result 
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of the leading order calculation obtained, by setting 77 = 0, in which case our 
equations reduce to the case studied inE3. 

We expect errors in our calculation to be of the order (ro/a)^, (kro)^ com- 
pared to the leading order. These errors are smaller than the experimental 
uncertainty in the finite energy case and of the same order as the experimental 
uncertainty in the case of the more precise measurement nearJe = 0, where we 
find ^ath = 6.33 ± 0.10 fmS compared to = 6.35 ± 0.02 fmi. 




Figure 3: k cot<5 in the J = 3/2 channel to ordepi (ro/a)" (dashed line) arfd (ro/a)^ (solid 
line). Circles are from the phase shift analysis inci and the triangle is fromEI. 

Our results seem to deviate from a simple effective range type expansion 
only around the pole at ^ 0.05 fm^'^ . (A pole in k cotS corresponds to a zero 
in the scattering matrix, which does not carry any special meaning.) This pole 
does not appear in potential model calculations (e.g.^)), and presumably will 
be smoothed out by higher-order terms that we have not yet included. It is 
interesting that the only "experimental" point in this region seems to indicate 
some structure there. More experimental information on this region would 
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be highly desirable to confirm this prediction. More experimental points are 
available only beyond break up and are not shown in Fig. ^. Notice that 
the break up region is still well within the range of validity of the pionless 
effective theory. It would be interesting to check how well the predictions 
would compared to experiment in this region. 

The calculation of higher-order corrections involves the knowledge of fur- 
ther counterterms like the ones giving rise to p-wave interactions, etc. These 
parameters can be determined either by fitting other experimental data or by 
matching with another effective theory — involving explicit pions — valid up 
to higher energies. If more precise experimental data — particularly at zero- 
energy — appear, we would be facing a unique situation where precision cal- 
culations in strong-interaction physics can be carried out and testedli3. There 
is the prospect of better measurements at zero energy actually being carried 
out by using neutron interferometry techniques, so this is not as far fetched a 
possibility as it may seem. 



5 The J = 1/2 channel 

The strong suppression of the three body force effects in the J = 3/2 channel 
is a reflection of the fact that, due to the Fermi principle, the three nucleons 
cannot occupy the same position in space. In the J = 1/2 channel, as well 
as in the case of three identical bosons, this restriction does not apply. One 
might expect that the J = 1/2 three body force term does not have to include 
two derivatives like the J = 3/2 three body force does and that, as a conse- 
quence, its effect on deuteron-nucleon scattering is suppressed by only (ro/a)'' 
(as opposed to (ro/a)^ in the J = 3/2 channel). There are however more 
profound differences between these channels. Let us consider, for simplicity, a 
system of three identical bosons such that there is a two boson bound state 
close to threshold. The leading order contribution is given by the graphs in 
Fig. H, with ro = 0. Power counting shows that all the graphs in this infinite 
series are ultraviolet finite. A closer look shows however that the sum of all 
graphs contains some ultraviolet divergence. This has been known for a long 
iime and diffeppt aspects of this phenomenon lead to the so called Thomas 
llJ and Efimovlla effects. Maybe the easier way to see this is to start from the 
analogue of equation (ffl) corresponding to the bound state of three bosons 



2 

tip,k) = - / dq t{q,k) 



pq^ 



ME\ 
ME J ' 



(20) 
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where E = -B + 3fc2/4M is the total energy. Equation (|2^) differs form 
only in the absence of the inhomogeneous term and in the coefficient in front 
of the kernel, that has a crucial opposite sign. We can map ( |2C| ) into a two 
dimensional quantum mechanics problem through the transformation 



X{r,p)^ I dpt(j>,k)^ : g-ry^-iu^^ (21) 

x(r, p) satisfies 



9;^ + Q^+MEjx{r,p)^0, (22) 
with the boundary conditions 

X(r,0) = (23) 

§^xir = 0,P)- Ixir = 0,p) + ^'xC^P, \p) = 0- (24) 

This same Schroedinger equation with these boundary equation would be ob- 
tained by considering the three particle problem in the usual wave mechanics 
formalism and using the Jacobi coordinates. The variable r corresponds to 
the distance between two of the particles and p to the distance between the 
third particle and the center of mass of the other two. Using polar coordi- 
nates R = yr^+p^, tga = r/p, we can see that this complicated boundary 
condition is simpler for R ^ a. Equation ( [2^ ) becomes 

^x(^,a = 0) + -^x(i?,V3) = (25) 

In the region i? <JC a we have then a separable problem with the solution 

X(i?, a) = V F,{R)sm{s,{TT/2 - a)) (26) 



where the Si are the solutions of 

g 

Si cos(si7r/2) = — ^sinsi7r/6. (27) 
v3 

Equation ( p7| ) has one imaginary solution isQ, with sq — 1.006. The equation 
for the radial part Fq{R) contains an attractive 1/R^ potential 
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Equation ( p8| ) is valid for 1/A ^ i? ^ a, where A is an ultraviolet cutoff. 
Taking a — > with A fixed we find that an infinite number of bound states will 
appear at threshold (Efimov effect). Keeping a fixed and taking A ^ oo bound 
states appear with arbitrarily low energy. That suggests that even at leading 
order there is need for a counterterm at r = p = (three body force). What 
is far from clear is whether only one counterterm would be enough to absorb 
the divergence or an infinite number of them (a whole form factor) would be 
necessary even at leading orderllj. Of course this last possibility would destroy 
the predictive power of the effective theory program in nuclear physics, at least 
on schemes in which pion exchanges are perturbative. 

Notice that the fact that three body forces appear at leading order in our 
approach does not contradict the usual statement that in phenomenological 
potential models they are small. By changing the cutoff the effect of two body 
forces can be transfered to the three body forces (in other words, the beta 
function for the three body force depends on the two body force). 

Another way of seeing that the ultraviolet problem described here cannot 
be seen diagram by diagram is to multiply the kernel of ( pO| ) by a parameter 
A. Perturbation theory in A at any finite order corresponds to the truncation 
of the series in Fig. ||. Repeating the arguments of this section we now find 
that the value of sq depends on a non analytic way on A. Sq is zero for A 
smaller than some critical value, and becomes non zero for larger values of A. 
Perturbation theory around A = will see any no hint of a finite sq and no 
attractive 1/i?^ potential. 
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